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ABSTRACT 

We propose a new approach to test possible corrections to Newtonian gravity using solar physics. 
The high accuracy of current solar models and new precise observations allow to constrain corrections 
to standard gravity at unprecedented levels. Our case study is Eddington-inspired gravity, an attrac- 
tive modified theory of gravity which results in non-singular cosmology and collapse. The theory is 
equivalent to standard gravity in vacuum, but it sensibly differs from it within matter, for instance 
it affects the evolution and the equilibrium structure of the Sun, giving different core temperature 
profiles, deviations in the observed acoustic modes and in solar neutrino fluxes. Comparing the pre- 
dictions from a modified solar model with observations, we constrain the coupling parameter of the 
theory, \k 9 \ < 3 ■ 10 5 m 5 s _2 /kg. Our results show that the Sun can be used to efficiently constrain 
alternative theories of gravity. 



1. INTRODUCTION 

In the last century General Relativity passed several 
stringent tests and it is now accepted as the standard 
theory of gravity and one of mankind's greatest achieve- 
ments (|Willll2005[ ). In the weak-field regime, General Rel- 
ativity reduces to Newtonian gravity, which is encoded 
in the famous Poisson equation for the gravitational field 



V 2 $ = AirGp, 



(1) 



where G is the gravitational constant and p is the matter 
density. In vacuum, the gravitational field of a spheri- 
cally symmetric mass M simply reads 



$(r) = -GM/r. 



(2) 



The validity region of the equation above has been tes ted 
and confirmed from submil limeter dHovle et al.ll2001l ) to 
solar system experiments (jWilll [2005l v! However, much 
less is known about Poisson's equation ((T|) inside mat- 
ter. In fact, the coupling to matter is one of the most 
delicate points in Einstein's theory. Several alternative 
theories have been proposed, which introduce modifica- 
tions in the coupling between m a tter a nd gravity (see 
e.g. iDamour fe Esposito-Faresel (|1993| )). The investi- 
gation of possible alternatives to the General Relativity 
paradigm are important. Extrapolating Einstein's the- 
ory to regimes in which it is not well-tested may lead to 
bias, potentially affecting astrophysical observations and 
our understanding of the Universe. 

At relativistic level, corrections in the gravity-matter 
coupling would affect the interior of neut ron stars and 
the c osmological evolution of the universe (jClifton et al.1 
120 111) . However, the uncertainty on the correct equation 
of state (EOS) describing the interior of a neutron star 
makes it difficult to disentangle the effects of an alterna- 
tive theory from those due to a different EOS. 



On the other hand, deviations from standar d gravity 
have been proposed even at N ewtonian level ([Milgroml 
[19831 iBanados fc Ferreir^[20l0l ) in a way which is com- 
patible with current experimental bounds. Theories such 
as these are consistent with all observations and at the 
same time are able to avoid long-standing problems of 
standard gravity. Thus, modified theories should be 
taken seriously and as important alternatives to explain 
our Universe and it is of utmost importance to develop 
methods to test and constrain them against standard 
gravity. 

In this work we propose a new approach, which is 
not affected by the degeneracy problems in neutron star 
physics and is complementary to cosmological tests. We 
shall investigative how deviations in Eq. (fTJ would af- 
fect the evolution and the equilibrium structure of the 
Sun and other stars, leaving potentially observable ef- 
fects. The high accuracy obtained with current stan- 
dard solar models and precise observations of the acous- 
tic modes and neutrino fluxes allow to perform stringent 
tests of the physics governing the sta r evolution and inte- 
rior (jTurck-Chieze fc Couvid at 2011). In the past, stellar 
evolution has been used to constrai n a possible t ime de- 
pendence of Newton's constant G ([Teller! 119481 ) . More 
recently, similar id eas have been used to put constraints 
on the value of G (lLopes fe Silk 2003 1), on the properties 
of dark matter particles (lLopes et alj|2002t lLopes fc Silkl 



l2010t iCasanellas fc Lop"esl 201 ID. and on the couplings of 
other particles ((G^ndolo_&Raffelt 2009). Finally, possi- 
ble modifications to the stellar st ructure in some alter- 
native scena rios were studied by iBertolami fc Paramos! 
((20051 [20081) using polytropic models. Given the high 
(and increasing) accuracy of present realistic solar mod- 
els and related observations, using the Sun as a theoret- 
ical laboratory where alternative theories of gravity can 
be challenged, is a very promising tool to constrain de- 
viations from Newtonian gravity. 
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ilidio.lopes@ist.utl.pt, vitor.cardoso@ist.utl.pt 
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2. PARAMETRIZED POST POISSONIAN 
APPROACH FOR MODIFIED GRAVITY 

The Parametrized Post Newtonian approach proved to 
be extremely efficient to constrain weak-fiel d deviation s 
from General Relativity in orbital motion (| Willi 120051 ) . 
The approach is based on a very general paramctriza- 
tion of the metric functions, and does not require any 
knowledge of the underlying alternative (metric) theory. 
Following a similar approach, here we parametrize vi- 
able couplings between matter and gravity in the non- 
relativistic limit, i.e. within Newtonian theory. We re- 
quire a modified Poisson equation which reduces to the 
usual one in vacuum, but which can accommodate extra 
terms in the coupling with matter. Assuming this theory 
is the non-relativistic limit of some covariant relativistic 
theory, we also require spacial covariance. Finally, we 
assume the theory contains at most second-order deriva- 
tives in the fields, although this condition can be easily 
relaxed. A general modified Poisson equation, up to sec- 
ond order in $, p and derivatives, which satisfies these 
requirements, reads 

V 2 $ = iwGp + *y v2 P + a g e ij Wi^VjP 
+np 2 + 7V/5 • Vp + eiV$ • Vp 
+e 2 $V 2 p + e 3 pV 2 $ + ... (3) 

The first term on the right hand side of the equation 
above is the standard Poisson term. The second one, 
proportional to K g , arises from the Eddington-inspired 
ravity theory recently proposed by iBanados &: Ferreiral 
20101 ). The other terms are higher order corrections and 
e y is the Levi-Civita symbol. All the parametrized cor- 
rections vanish in vacuum, so that the theory above is 
consistent with the inverse square law behavior (|2"j), but 
most of the extra terms in Eq. ([3]) violate the equivalence 
principle and ar e therefore already strongly constrained 
by experiments (|Willll2005() . Two notable exceptions are 
the terms proportional to n g and, for spherically symmet- 
ric configurations, the term proportional to a g . These 
two terms are consistent with the equivalence principle, 
and mostly unconstrained presently. 

A CASE STUDY 

For concreteness, here we focus on a particular case, 
setting 7 = rji — e, = in Eq. fl3J|. The modified Poisson 
equation reduces to 



V 2 $ = AwGp + -j-V 2 p + a g e ij Vi<S>V 3 p, 



(4) 



where 



] = cm 5 /(g s 2 ) = [[G]][[i? 2 ]]. Requiring spher- 
ical symmetry, the hydrostatic equilibrium equation fol- 
lows 

dP Gm(r)p 



dr 



-pp . 



(5) 



where no terms proportional to a g arise, due to the spher- 
ical symmetry. The choice 7 = t]i = e = is moti- 
vated by several reasons. First of all, the terms we are 
neglecting would introduce violations to the equivalence 
principle, which is e xperimentally confirmed with great 
precision (Will 2005|). Secondly, the equation above rep- 
resents the most general modified Poisson equation which 



is first order in $, p and satisfies the requirements pre- 
viously discussed. The extra terms would only intro- 
duce higher order corrections. Furthermore, this the- 
ory is the nonrelativistic limit of a well-motivated theory 
of gravity, which prevents the formation of singularities 
in cos mology and in the stellar collapse of comp act ob- 
jects (jBanados fc Ferreiral [20ToHPani et al.ll2011[) . Here 
we investigate how this theory would modify the interior 
and the evolution of the Sun. 

3. THE EVOLUTION OF THE SUN 

The high accuracy of current solar models and precise 
observations allow to test standard gravity against al- 
ternative theories at unprecedented levels. Alternative 
theories would affect the evolution and the equilibrium 
structure the Sun, giving different core temperature pro- 
files and deviations in the observed acoustic modes and in 
solar neutrino fluxes. Comparing the predictions from a 
modified solar model with observations, we can constrain 
the coupling parameter of the theory. 

Modeling the solar interior not only requires to de- 
scribe the present solar structure, but also to ex- 
plain the evolution of the Sun from the ignition 
of hydrogen nuclear fusion to the present day (sec 
e.g. iTurck-Chieze fc Lopesl (119931 ) for a review). The so- 
lar models are constructed on the basis of plausible as- 
sumptions, which translate in a set of four ordinary dif- 
ferential equations. The star is considered in hydrostatic 
equilibrium, which means that the hydrostatic pressure 
resulting from the thermonuclear fusion of hydrogen to 
helium must be exactly balanced by gravity. The nuclear 
reactions are produced in the pp chain and in the CNO 
cycle, the former affecting more strongly the tempera- 
ture profile in the solar core. Furthermore, the star is 
assumed to be in thermal equilibrium, i.e. the energy 
produced by nuclear reactions balances the total energy 
loss via radiative energy flux and via the energy carried 
away by neutrinos. 

Within ~ 70% of the solar radius, the most efficient 
transport mechanism of energy from the solar center out- 
wards to the stellar surface is due to electromagnetic ra- 
diation, while in the outer region, the so-called convec- 
tive zone, the energy is mainly transported by convec- 
tion. The radiative energy transport (and, in turn, the 
temperature profile) is governed by the Rosseland mean 
opacity, which takes into account that photons interact 
with electrons and ions in the dense plasma in the so- 
lar interior, while they mostly interact with atoms and 
molecules at the solar surface, where radiative transport 
is again significant. 

One of the basic assumptions of any solar model, 
namely the hydrostatic equilibrium, ultimately depends 
on how strong and efficient the gravitational self- 
interaction is inside the Sun, i.e. on Poisson's equa- 
tion ([3]). Any corrections would affect the thermal bal- 
ance and, in turn, the temperature profile inside the star, 
leaving potentially observable signat ures. 

Fin ally, effects due to rotation (|Pinsonneault et al.l 
119891 ) and magnetic fields (jPassos fc Lopesll2008l ) are usu- 
ally neglected in standard solar models. These processes 
take place on a much shorter timescale than the evo- 
lutionary timcscalc of the Sun and their inclusion re- 
sult s in minor structure chang es in the solar interior (see 
e.g. ITurck-Chieze et al.l pOTO )). 
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3.1. Equations governing stellar equilibrium and 
evolution 

Under the previous assumptions, the internal structure 
of the Sun is governed by the following ordinary differ- 
ential equations for (r, P, L,T), the radius, pressure, lu- 
minosity and temperature, respectively 



dr 
dq 
dP 
dq 
dL 
dq 



Airr 2 p ' 
GM%q 
47rr 4 



■■Mr, 



k„ dp 
T P d~q 



dS_ 
~dt 



(6) 
(7) 
(8) 



where q — to/A/q is a convenient choice of the indepen - 
dent variable, since mass loss is ncg lected (|Clavtonlll968D . 
The first and third equations above are the standard con- 
tinuity equation and conservation of thermal energy, re- 
spectively, whilst the second equation describes the hy- 
drostatic equilibrium (note that Eq. ([7]) is equivalent to 
Eq. ([5]) when expressed in terms of independent variable 
q). Finally, the equations above must be supplied by an 
appropriate transport energy eq uation, for th e convective 
zone and for the radiative zone ([Morel 1993). Due to the 
modified Poisson equation ([7]) , the standard equation for 
the convective energy transport is indirectly modified as 
follows 



dT _ dP dT 
dq dq dP 



GM%q 



inr 4 



dp 
t P dq-. 



T 



(9) 



where V = (ilogT/dlogP is the temperature gradient. 
For adiabatic changes, the temperature gradient can be 
simply relate d to one of the adi abatic exponents, V a d = 
(r 2 - l)/r 2 (jWeiss et all 12004) . In the radiative zone, 
the transport energy equation is unaffected by n g and it 
simply reads 



dT_ 
dq 



3M k L 



(10) 



where k is the Rosscland mean opacity and a is the Boltz- 
mann constant. 

3.2. Numerical procedure 

The modified equations above have been included, to- 
gether with all_the relevant physical processes, in CE- 
SS AM (Morel 1997), a self-consistent numerical code for 
stellar structure and evolution. The main physical in- 
puts of the solar models are the fo llowing: the nu- 
clear reaction rat es are t aken from [Adclbcrgc r et all 
(|1998f ). with the iMitlerl (|1977| ) intermediate screen- 
ing; the opacities are taken from the OPAL95 ta- 
bles (jlglesias fc Rogers] 11996ft for tempera t ures a bove 
5600 K and from I Alexander fc Ferguson! (|1994D for 
lower temperatures; we used the tabulated OPAL 
EOS (jRogers et all I1996D ; microsc opic diffusion is in- 
cluded following the prescription of TMichaud fc Proffittl 
(Il993f): finally, the s olar abundances are taken 
fromlAs plund et all (|2005l ). For comparison, all the mod- 
els were also com puted with the older, low -metallicity 
solar abundances ([Grevesse fc Sauval I1998D . leading to 
nearly identical results. This test stresses the robustness 
of our approach and shows that our analysis is virtually 



independent of the current uncertainties of solar model- 
ing. 

In order to constrain the values of the coupling pa- 
rameter Kg which are compatible with present observa- 
tions of the Sun, we constructed calibrated solar models 
for different values of n g . The models are calibrated to 
fit the solar properties with an accuracy of 10~ 5 . The 
calibration is performed by varying the parameters Xq 
(the initial abundance of hydrogen in the young Sun) 
and a (which parametrizes the efficiency of convection 
as a mechanism of energy transport), and by fixing the 
solar surface heavy-element content (Z/X)q = 0.0165, 
age t Q = 4.57 Gyr, radius R = 6.9599 x 10 10 cm, 
mass M Q = 1.9891 x 10 33 g and luminosity L Q = 
3.846 x 10 33 ergs" 1 . 

It was possible to construct calibrated solar models for 
-0.032Gi?| < K g < 0.02Gi?|. For k 9 < -0.032Gi?|, 
no equilibrium st ars can be constru cted, in agreement to 
what is shown in IPani et all (|2011D for simple polytropic 
models. On the other hand, for K g > O.Q2GR Q equilib- 
rium stellar configurations can be constructed, but their 
internal structure is so strongly modified that the ob- 
served solar properties (namely Z/X, t@, i? Q , Mq and 
Lq) cannot be matched simultaneously. In Table [TJ we 
show the values of Xq and a required to calibrate the so- 
lar models, together with the central temperature, den- 
sity and pressure of the models. 



Kg 


X 


a 


T 


Pc 


Pc 


(GR%) 






( lO^K ) 


(g cm -3 ) 


( dyn cm -2 ) 


-0.032 


0.78 


2.77 


15.54 


161.5 


2.56X10 17 


-0.01 


0.76 


2.07 


15.25 


150.1 


2.36X10 17 





0.74 


1.84 


15.18 


146.7 


2.29X10 17 


0.01 


0.73 


1.65 


15.12 


144.8 


2.25X10 17 


0.02 


0.72 


1.48 


15.09 


143.6 


2.22X10 17 



TABLE 1 

Characteristics of some solar models for different values 
of Kg. All stellar models have M = Mq, L = Lq, and 

R = Rq AT THE SOLAR AGE t n = 4.57 GYR. 



Solar models with K g > have a lower central density 
and a lower core temperature, whereas models with K g < 
work in the opposite direction. These results can be 
qualitatively understood as follows. Equation ([S]) can be 
written in a more evocative form as 



dP m(r)p(r) 
— = -G cS (r)— — 



(11) 



where we have defined an "effective" Newton's constant 

K g r 2 p'{r) 



Goff(r) = G 



4 m(r) 



(12) 



Since p'{r) < inside the Sun, G e g ^ G when n g ^ 0. 
When K g < 0, we expect a stronger effective gravitational 
force which, for main sequence stars in hydrostatic equi- 
librium, leads to an increase in the central temperature 
and, consequently, in the rate of thermonuclear reactions. 
The solar models constructed with K g < also require 
a higher initial abundance of hydrogen to match the so- 
lar observables. This fact can be explained by homology 
scaling (see e.g. iKippenhahn fc Weigertl (|1990ft ): the lu- 
minosity of a star scales as a high power of G and the 
mean molecular weight, therefore an effective increase in 
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Fig. 1. — Density profiles of the modified solar models computed Fig. 2. — 8 B neutrino flux predicted by our modified solar models 

with different values of K g . normalized to the flux predicted by our standard solar model. 



the gravitational force must be compensated by a de- 
crease in the mean molecular weight (and consequently 
an increase in the hydrogen abundance) to achieve the 
same luminosity. This qualitative picture is in agreement 
with results ob tained for diff e rent ( constant) values of G 
(cf. Table 1 in lLopes fe Sitkl (l200l ). 

4. RESULTS 

The Eddington-inspired theory of gravity leads to 
strong modifications in the solar structure. In a wide 
region of the parameter space of the theory, the modi- 
fied solar models show important variations in the cen- 
tral temperature and in the density profile (see Fig. [T]). 
These two signatures can be tested against solar observ- 
ables. In particular, we shall show that solar neutrino 
measurements (which are sensible to T c ) and helioseis- 
mic acoustic data (sensible to sound speed and density 
profiles) strongly constrain the values of the parameter 
Kg that are compatible with present observations. 

4.1. Solar neutrinos 

Solar neutrinos provide a unique window to the so- 
lar interior due to the high sensitivity of thermonu- 
clear reactions to the temperature at which they 
take place. In particular, the 8 B flux, produced in 
the inner 10% (in radius) through the pp chain, is 
very sensitiv e to the central temperature o f the Sun: 
08 B oc T c 18 (ITurck-Chieze fc Couvidatl l20ll . The pre- 
dicted neutrino flux is expected to depend strongly on 
Kg since different couplings lead to different central tem- 
peratures (modified solar models may lead to variations 
of up to 3% in T c ). The observed 8 B neutrino flux is 
currently measured with high precision by neutrino tele- 
scopes: (5.046 ± . 16) x 10 6 cm~ 2 s -1 (|Aharmim et al.l 
I201tt iBellini et al.ll2010l ). Thus, the theory can be con- 
strained on the basis of incompatibility with observa- 
tions. 

Our results for the solar neutrino fluxes are shown in 
Fig. [5J As expected, the dependence on the coupling- 
parameter Kg can be understood in terms of effective 
gravitational constant. Positive values of K g lead to a 
smaller G e s, contributing to a lower central temperature 
and, in turn, to a lower expected neutrino flux. Negative 
values of K g work in the opposite direction. 

The theoretical uncertainty of standard solar model- 
ing have to be consistently taken into account when 
comparing the predictions of our models with the ob- 
servations. Previous works have shown that the largest 



source of uncertainty in the calculation of the solar 
neutrino fluxes comes from the uncertainty in the val- 
ues o f the surface heavy element abundances of the 
Sun dBahcall et al l 120061: iGonzalez-Gartia et al.l 120101: 
INorena et al.l 120 111) . iBahcall fc Serenellil (120051) deter- 
mined, using Monte Carlo simulations for 10000 solar 
models, that the total ler theoretical uncertainty in the 
predicted 8 B neutrino flux is below 17%, in the most con- 
servative scenario. In addition, we also take into account 
the deviation of 20% in the predicted 8 B flux when differ- 
ent esti mations of the so l ar abu ndances are implemented 
(see e.g. iSerenelli et al.l ([2009D ). Considering both the 
theoretical and experimental uncertainties, we estimated 
that models that predict a 8 B flux which deviates more 
than 30% from our standard solar model can be conserva- 
tively ruled out, in agr eement with the t hreshold consid- 
ered by other authors ([Taoso et al.H20l6T) . Following this 
analysis, we conclude that values of K g < 



0.024Gi?| 



are excluded by the observation of 8 B solar neutrinos 
(see Figure [2]). 

On the other hand, we found that the 7 Be neutrino 
flux only provides less stringent constraints on K g , as 7 Be 
neutrinos are produced in a wider region in the center of 
the Sun and, consequently, a re less sensitive to its cen- 
tral t emperature (4> 7 Bc oc T 8 (jTurck-Chieze fc CouvidatJ 

Enni)). 

4.2. Helio seismology 

The solar acoustic modes are nowadays measured 
with exquisite precision by helioseis mic missions aboard 
space crafts, such as GOL F/SOHO (fTurck-Chieze et ahl 
E993), M DI/SOHO (IScherrer et al.l 119951 ) and 
HMI/SDO dZhao et al.l 120 111), and by ground net- 
works such as BiSON jBroomhall et al.l l2009h and 



GONG ([Harvey et al.|[l996 ). The analysis of helioseismic 
data has provided a valuable tool to probe the solar inte- 
rior, revealing the sound-sp eed and density profiles down 
to 10% of the solar radi us (jChristensen-D alsgaar d et al.1 
[TMllGough et al.iri996h . 

Different helioseismic parameters have been 
used to investigate v a rious aspects of sola r 
physics ([Thompson et al.l 119961 : iGizon et al.l 120101 ). 
In particular, the small separation between the frequen- 
cies of modes with different degree I and radial order 
n, 5v n j = v n j — v n -i,i+2i is a helioseismic quantity 
which is very sensitive to the te mperature gradient 
in th e deep interior of the Sun (lOti Floranes et all 
120051) . In addition, the modes with degree I = 



■5 




1.5 2 2.5 3 

v ( mHz ) 




-0.03 -0.02 -0.01 0.01 0.02 



K S ( G R|) 

Fig. 3. — a) Small separations for 1=0 calculated in o ur models 
compared with helioscismic data (Broomhall ct al. 2009). b) Mean 
small separation for 1=0 and v > 2000/^Hz for our modified solar 
models, normalized to the prediction for re g = 0. 

correspond to acoustic waves that traveled through 
the entire stellar radius and ca rry information about 
the density profile of the Sun (ILopes fe Turck-Chiezd 
[1991 [RmdTurgh fc Vorontsovl 12000ft . Therefore. 8v„ ,-n. 
also known as fine spacing or c?02, is a very suitable 
parameter to detect the signatures that alternative 
theories of gravity leave on the solar interior. 

The small separations in modified solar models are 
compared with solar data in Figurc[3]a). As expected, for 
Kg = the fine spacings exhibit a moderate disagreement 
with the observations. This discrepancy, which disap- 
pears when the older, low-Z solar abu ndances are consid- 
ered, has been discussed in detail bv iBasu et al.1 <|2007f ). 
On the other hand, for large values of K g the deviations 
from helioseismic data are much larger, providing a clear 
way to discriminate viable models. 

As discussed for solar neutrinos, when the helioscismic 
quantities are used to constrain solar models in modified 
theories of gravity, the uncertainties of solar modeling 
also have to be taken into account. Compared to so- 
lar neutrinos, the theoretical uncertainties on the mean 
small separation {5v n ,i=o) are much smaller. The varia- 
tion on (5i>n,i=o) res ulting from differe nt solar models is 
of the order of 2-3% (|Basu et al.ll2007t) . Considering this 
uncertainty, we can rule out those models that lead to 
deviations in (Sv n ,i=o) greater than 4%. This diagnos- 
tic establishes strong constraints on n gi ruling out the 
regions n g > 0.016 G R| and n g < -0.01 G R|. (sec 
Figure Ob). 

4.2.1. Other helioseismic constraints 

Another constraint on deviations from Newtonian 
gravity comes from the comparison of the solar and 
model sound speed profiles, the former being ob- 
tained with high precision from helioseismic observa- 
tions. Remarkably, the standard solar model repro- 
duces the sound speed profile of the Sun with an ac- 
curacy better than 1% in most of its interior. How- 
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r(R G ) 




-0.03 -0.02 -0.01 0.01 0.02 



K g ( G R|) 

Fig. 4. — a) Relative differences between the sound speed pro- 
files of our modifi ed solar models and th e solar sound speed from 
helioseismic data (Broomhall et al. 2009). b) Mean deviation be- 
tween the solar and model sound speed profiles, normalized to our 
standard solar model. 

ever, right below the convective envelope the devia- 
tions from the observed sound speed are larger (this 
discrepancy is common in mo dels adopting the lat- 
est, high-Z solar abunda nces dMontalban et al] 120041 : 
iDelahave fc PinsonneauTtl 120061: iSerenelli et all 12009ft ). 
Figure |4]a) shows the relative differences between the 
helioscismically inverted and the sound speed profiles 
of some of the modified solar models, dc/c = (c© — 
c-modei) I Cmodei- The mean difference (\Sc/c\) is a measure 
of how accurately a solar model reproduces the sound 
speed profile of the Sun, and therefore it can be used to 
put constraints to modified theories of gravity. Those 
models leading to a relative deviation (|<5c/c|) more than 
two times larger than the (\Sc/c\) obtained for Kg = can 
be conservatively ruled out. As shown in FigurelUb), the 
constraints from the sound speed profile rule out models 
with K g > 0.012 G Rq. 

Hclioscismology also provides accurate measurements 
of the depth of the convec tive envelope, Rcz = 0.713 ± 
0.001 (jBasu fc AntiallT997T ) and the helium surfa ce abun- 
dance Y s = 0.2485±0.0035 (jBasu fc Antial[200l . Monte 
Carlo simulations have shown that the theoretical un- 
certainty fr om solar modeling i s below 2% for Rcz and 
5% for Y s (jBahcall et all [2001 . Consequently, we can 
conservatively rule out models that predict deviations of 
these quantities larger than 3% and 7%, respectively . 
As shown in Figure EJ this allows us to put the following 
constraints on the parameter K g : — 0.016 GRq < n g < 
0.013 GR% and n g > -0.018 GR%, respectively from the 
observations of Rcz and of Yg- 

5. DISCUSSION AND CONCLUDING REMARKS 

Our results show that the Sun is a very good testing- 
ground to constrain generic modified theories of gravity, 
for instance theories such as the ones described in Eq. §§§ 
and even more exotic or yet to be proposed corrections. 
For the particular case of Eddington-like theories, Tablc[2] 
summarizes the constraints on the coupling parameter n g 
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Fig. 5. — a) Depth of the convective envelope Rcz and b) helium 
surface abundance Y$ of the modified solar models normalized to 
our standard solar model. 



of the theory. Our results show that, in order to obtain 
a viable solar model, a theory as general as Eq. Q is 
strongly constrained. Combining all the experimental 
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Range of K g excluded 
-0.024 G R| 



< 



Kg < 



-0.01 G R| and K g > 0.016 G R| 
K g > 0.012 G R| 
< -0.013 G R2 and k„ > 0.013 G R%, 





<-0.018 G R; 



G 



TABLE 2 



USING DIFFERENT SOLAR CHARACTERISTICS. 



bounds, the coupling constant n g must lie in the interval 
-0.01 < k 9 /(GRq) < 0.012, i.e. approximately |k b | < 



10 5 m 5 E 



7kg 



It is important to stress that this result does not rule 
out Eddington-inspircd theory as a promising alterna- 
tive to Einstein's theory. Previous studies showed that 
most of the appealing features of the theory would persist 
eve n for a (positive) arbitrarily small coupling parame- 
ter (jBanados fc Ferreiral |20T1 IPani et all [20111) . which 
is perfectly consistent with current observations of solar 
neutrinos and helioseismology. 

Modified gravity is also relevant as an alternative ap- 
proach to the solar abundance problem. The particu- 
lar theory we considered only offers a partial solution 
to this problem. Indeed, models with K g < predict 
the base of the convective envelope at a smaller ra- 
dius than the standard solar model, reconciling the pre- 
diction with the hclioseismically inferred value. How- 
ever, the predicted helium surface abundance for the 
same models with n g < is then even more un- 
derestimated than for standard solar models. Simi- 
lar pa rtial solutions wer e discussed in different con- 
texts dCastro et all l2007t IChristensen-Dalsgaard et all 
20091: iGuzik fc Mussackll2010t ISerenelh et al. | |2011[) . Al- 
though Eddington-inspircd gravity suffers from the same 
limitations, other gravitational corrections could affect 
the solar interior in a different way and they should be in- 
vestigated more carefully. We leave this interesting topic 
for future work. 
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